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Abstract 

Thermal instabilities can cause a radiative shock to oscillate, thereby modulating the 
emission from the post-shock region. The mode frequencies are approximately quantised in 
\ analogy to those of a vibrating pipe. The stability properties depend on the cooling processes, 

the electron-ion energy exchange and the boundary conditions. This paper considers the 
effects of the lower boundary condition on the post-shock flow, both ideally and for some 



(N 



■^j- , specific physical models. Specific cases include constant perturbed velocity, pressure, density, 

flow rate, or temperature at the lower boundary, and the situation with nonzero stationary 
■ flow velocity at the lower boundary. It is found that for cases with zero terminal stationary 

velocity, the stability properties are insensitive to the perturbed hydrodynamic variables at 
the lower boundary. The luminosity responses are generally dependent on the lower boundary 
Qh' condition. 

o 
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'■ 1 Introduction 



A system that depends on the length- or time-scales of radiative cooling or other energy exchange 
processes may be thermally unstable (Field 1965) if perturbations from the initial or equilibrium 
condition alter the cooling scale in a manner that enhances the perturbation. The effects of 
radiative cooling are expressed in terms of a function of local hydrodynamic variables such as 
density and temperature. Qualitatively, a thermally unstable system has a form of cooling function 
which increases the cooling length when the shock temperature increases and conversely the cooling 
length decreases when the shock temperature lowers. Langer, Chanmugam & Shaviv (1981) 
discovered such a thermal instability in the radiative accretion shocks of white dwarfs (see e.g. 
reviews by Cropper 1990, Wu 2000) This instability was verified by subsequent studies using 
different numerical techniques (e.g. Imamura, Wolff & Durisen 1984), and linear stability analyses 
(e.g. Chevalier & Imamura 1982). 

Linear analyses (Chevalier & Imamura 1982; Imamura et al. 1996; Saxton et al. 1998; and 
Saxton & Wu 1999) indicate that the shock oscillations have a sequence of eigenmodes. If 
bremsstrahlung dominates the cooling then there tends to be a fundamental mode stable against 
oscillations, and unstable overtones. Depending on its relative efficiency, cyclotron cooling tends 
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to stabilise modes, except when the electron-ion energy exchange is comparatively inefficient, in 
which case modes can be destabilised with increasing cyclotron efficiency. 

Toth & Draine (1993) found that when bremsstrahlung cooling is the dominant cooling process 
then the frequencies resemble those of a pipe open at one end, (oscillation frequency / oc n — | for 
integer harmonic number n). This pattern persists for two-temperature shocks when the electron- 
ion exchange is efficient (Imamura et al. 1996, Saxton & Wu 1999, Saxton 1999) When cyclotron 
cooling is very efficient, so that the electron- ion exchange process is unable to maintain strong 
thermal coupling between electrons and ions, the modes are more like those of a doubly-closed or 
doubly-open pipe / oc n, (Saxton & Wu 1999). 

The standard analytic treatment defines the lower boundary as the place where the flow ve- 
locity becomes zero in the time-independent solution and does not oscillate in the time-dependent 
response to perturbation (e.g. accreted material diffuses across a stellar surface very slowly com- 
pared to the oscillatory timescale). The pressure, density and other hydrodynamic variables are 
nonzero and are allowed to oscillate. Different numerical investigations treated the lower boundary 
differently: e.g. dense matter near the lower boundary was deleted in Imamura (1985); was merged 
with neighbouring cells in Imamura, Wolff & Durisen (1984); and may have had its gradients soft- 
ened by the spatial filter on the Eulerian grid of Langer, Chanmugam & Shaviv (1981). These 
diverse treatments produced some different results. Langar et al. (1981, 1982) find a fundamental 
mode that is unstable to oscillations; whereas other studies (e.g. Imamura 1985; Imamura et al. 
1984; Wolff et al. 1991; Wood et al. 1992) found that the first unstable mode is the first overtone, 
and the fundamental mode is damped, consistently with stability analyses. The collaborations also 
produce different results regarding the phasing of oscillatory modulations of the bremsstrahlung 
and cyclotron luminosities. 

This paper investigates the role of the lower boundary in relation to the stability and emission 
properties of post-shock flows. This issue has not been thoroughly discussed in previous studies of 
one-temperature (Saxton et al. 1998) and two-temperature (Saxton & Wu 1999, 2001) radiative 
shocks. In an appendix, Imamura et al. (1996) stated that changing the lower boudnary condition 
from constant-velocity to constant-pressure did not affect their results with a power-law cooling 
function. The present work explores different analytic choices of the lower boundary condition, 
in the presence of two cooling processes and inequality between the electron and ion tempera- 
tures. White dwarf accretion shocks are an illustrative case study. The conditions considered 
herein include constant pressure, density, flow-rate and temperature at the white-dwarf surface, 
or a boundary with nonzero longitudinal velocity (e.g. due to lateral "leakage" transverse to the 
magnetic field lines). 

The plan of the paper is as follows. Section |2| outlines the general formulation, and section [3] 
describes the cases of different boundary conditions investigated. Section |] presents and discusses 
the local and global stability properties. Section [| examines how the different lower boundary 
conditions affect the modulation of the luminosities of the oscillating post-shock flow. Section ^| 
draws the conclusions. 



2 perturbative analysis 

The hydrodynamics of the accretion flow are governed by the equations for continuity of mass, 
momentum, total energy and electron energy: 

JUvVU + P(V-v) = 0. (1) 
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where v, p, P and P c are the fluid velocity, density, total pressure and electron pressure, and 7 is 
the adiabatic index (assumed equal to 5/3 for an ideal monatomic gas). The radiative loss term 
A = A br + A cy comprises a sum of energy losses due to bremsstrahlung and cyclotron cooling. The 
functional form of A was described in Wu (1994). 



A = A br +A 2 
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where p s and P e ,s are the density and electron partial pressure at the shock, and the bremsstrahlung 
cooling function is A br = Ap 2 (P c / p) 1 ^ 2 with A m 3.9 x 10 16 in c.g.s. units (see e.g. Rybicki & 
Lightman 1979). The parameter e s is the efficiency of the second cooling process compared to 
bremsstrahlung cooling, evaluated at the shock. For the case of cyclotron-emitting post-shock 
regions of accreting white dwarfs, the indices are a = 2.0, (5 = 3.85. 

The term T is the electron-ion energy exchange, (e.g. Spitzer 1962, Melrose 1986, Imamura et 
al. 1996). 
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where m^ e , T^ e , 6- he = k B T^ e /mi te c 2 , and are the ion and electron masses, temperatures, 
dimensionless temperatures and number densities of ions and electrons respectively. The other 
constants are k B the Boltzmann constant, c the speed of light, e the electron charge, and InC is 
the Coulomb logarithm. 

Normalising the hydrodynamic quatities by their respective pre-shock values (vg the free-fall 
velocity, p a the pre-shock density of the flow, and x s o the equilibrium shock height), the variables 
of the stationary solution are defined as To = v/vg, no = P/p a Vg and 7r c = P e /p a Vg for the flow 
velocity, total pressure and electron pressure. The hydrodynamic variables are functions of the 
vertical position coordinate £ = x/x 8 , where x s is the instantaneous shock height, as distinct from 
the equilibrium value. Then the cooling function and electron-ion exchange function are: 



A 

r 



(p a Vs/x s0 )ilj c TQ /2 7T c 1/2 [1 + e s /(p a , P c )} 



(p a Uff/aJB0)^c^«T0 71 c 3/2 (1 - T - 2tT c 



(7) 
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where the constant ip c = Axgop^/vg and ip e \ parameterises the efficiency of electron- ion energy 
exchange compared to the radiative cooling. The conditions at the shock (£ = 1) are Tq — |, 
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where <r s = (P e /Pi) s is the electron to ion pressure ratio at the shock. 
The stationary state of the system is completely specified by the parameters (cr s , ^5 e s)- 

As in Imamura et al. (1996) and Saxton & Wu (1999), the shock position is perturbed about 



its equilibrium position, x s = x s0 + Xgie"', with complex frequency to = 5vg/x s0 . The stationary 
solution is obtained by separating the zeroth-order terms in the hydrodynamic equations, involving 
the total pressure tiq and electron pressure n c . 
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The perturbation of the post-shock structure is described by a set of differential equations, 
with the complex variables Af(£), A T (£), A 7r (^) and A e (£), corresponding to perturbations of den- 
sity, longitudinal velocity, total pressure and electron pressure respectively. The A functions are 
normalised by e = x s iu/vs, the relative amplitude of the oscillation of the shock position. The 
instantaneous density, velocity, total pressure and electron pressure are 

p{Z,y,t) = Pa-Co(0(i + ^A c (0^) 

v(£,y,t) = -^-To(£)(l + £A r (£)e"" 
P(£,y,t) = Pa v 2 s ■ n (0 (l + eXAOe^ 

Thus we obtain the first-order perturbed time-dependent equation 
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(see Saxton & Wu 1999), where A = (7 — l)A(x s0 / p a Vg) and r = —Vo/vg are dimensionless 
versions of the cooling function and accretion flow velocity. The Fi,...,Fs terms are complex 
functions involving the perturbed variables and the variables of the stationary solution. 

The differential equations are integrated from the shock down to the lower boundary. Values 
of 5 which satisfy the appropriate boundary conditions are the eigenfrequencies. The complex 
eigenvalues 5 = S R + i5 l describe the oscillatory properties of the shock in different modes. The real 
part, 5 R , is a growth term: if 5 R > then the shock is unstable to oscillations in the corresponding 
mode. The imaginary part, 5 1 , is proportional to the oscillatory frequency of the mode. Figure |l] 
shows examples of how the eigenplanes of the perturbed variables — log |A r | evaluated for a given 
set of boundary conditions. The A-profiles corresponding to the 6 eigenvalues are the eigenfunctions 
of the post-shock structure's local oscillatory response to the shock oscillations. 

For each cooling process, the luminosity response is calculated by multiplying the stationary- 
state cooling function and corresponding perturbation, and then integrating their product over 
the entire post-shock region. The luminosity response of a mode, is expressed in terms of the 
stationary-solution luminosity and normalised by the size of the shock-height oscillation, i.e. 
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where Ab r = §Af + |A e and A cy = — 2.35A^ + 2.5A e (see Saxton & Wu 2001). The absolute value 
of the ratio ( |T4"D is the relative amplitude of the bremsstrahlung(cyclotron) luminosity oscillation, 
and the phase $br (or $ cy ) is relative to the oscillatory phase of the shock motion. 



3 Boundary conditions 
3.1 shock jump condition 

The boundary conditions at the shock (£ = 1) are determined by the strong-shock Rankine- 
Hugoniot conditions. It can be shown that the stationary-solution conditions are tq = |, 7Tq = |. 
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Figure 1: Complex frequency planes showing the eigenmodes of systems with parameters 
(a s , ipei, e s ) = (0.5,0.5,0) in the top row and (er s , ip ci , e s ) = (0.5,0.5,100) in the bottom row. 
The lower boundary condition is A r = 0, tq = in the left column, and A r = 0, tq = 0.01 in the 
right column. The contour levels are at logarithmic intervals of the function 1/|A T |; the modes are 
at the peaks. The real part of the an eigenvalue, 5 R , describes the instability of a mode: positive 
indicates instability. The imaginary part of the eigenvalue, S I , is proportional to the frequency of 
oscillation. 
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If the ratio of electron and ion pressures is parameterised as a s then the electron pressure condition 
is n e = f/(l + o-s^ 1 )- The conditions on the perturbed variables are derived in the frame that is 
comoving with the shock (Imamura et al. 1996, see also appendices of Saxton et al. 1997), A^ = 0, 
A T = —3, and A^ = A o = 2. The system parameters are the same as defined in Saxton & Wu 
(1999): e s is the efficiency of cyclotron cooling relative to bremsstrahlung cooling at the shock; cr s is 
the ratio of electron to ion pressures at the shock; specifies the efficiency of electron-ion energy 
exchange relative to the total efficiency of radiative cooling. The shock boundary conditions are 
the same in all cases studied; alternative conditions at the lower boundary are considered. 

3.2 lower boundary conditions 

The post-shock flow has a lower boundary, labelled £ = 0, and the form and value of the boundary 
conditions there depend on the nature of the system under consideration. For example, the bound- 
ary of cooled gas far downstream from an interstellar medium shock and the boundary formed 
by a stellar surface in an accretion flow may impose different conditions on the hydrodynamic 
variables. 

3.2.1 "perfect" stationary- wall 

If the flow velocity vanishes at the lower boundary then the condition is r = at £ = 0. This is 
sometimes known as the "stationary-wall" condition. The additional choice of A r = means that 
the flow velocity is absolutely zero and does not oscillate at the lower boundary. The "perfect 
stationary-wall" boundary condition is the conventional choice in many stability analyses of radia- 
tive shocks (e.g. Chevalier & Imamura 1982, Toth & Draine 1993, Imamura et al. 1996, Saxton et 
al. 1998). Bertschinger (1986) demonstrated that this condition also applies to an ISM radiative 
shock, at the inner boundary which is defined as the place where gas is cooled to a point where 
total cooling is zero. 

Wolff, Gardner & Wood (1989) carried out numerical studies of accreetion onto white dwarfs 
where cyclotron cooling does not dominate, and bremsstrahlung cooling, Compton cooling, electron- 
ion energy exchanges and electron thermal conduction are included. They considered a "perfect 
stationary-wall" condition at the lower boundary, with additional constraints: (i) a constant 
nonzero electron temperature T e = T cut , the cut-off value where the total cooling A = 0; (ii) a zero 
gradient in temperature, dT e /dr = 0. Thes presence of conduction and the alternative form of 
A must provide a solution for time-independent flow structure that is significantly different from 
that considered in the present work, which focusses on the role of lower boundary conditions on 
the perturbed variables. The effects of boundary conditions on stability properties of structures 
where conduction and finite cooling cut-offs are physically important comprises a significant topic 
for future investigation. 

3.2.2 modified stationary-wall 

The A T = fixed-velocity assumption is not the only choice of lower boundary condition which 
is consistent with the stationary solution that provides r = at £ = 0. There are alternative 
physical propositions. 

The stationary-solution pressure at the lower boundary is a well-defined value, ttq = 1- The 
condition A^ = means that the total pressure of the flow does not oscillate at the lower boundary. 
This is a kind of hydrostatic condition. In the case of accreting white dwarfs, this condition would 
mean that the pressure wave of the shock oscillation does not enter the white dwarf atmosphere. 
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Table 1: The boundary conditions studied in this paper, in terms of the stationary and perturbed 
variables. 



case stationary variables perturbed variables note 
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To 


= o, 


7T = 
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A r = 


§P-2.1| "perfect stationary- wall" 
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§[3. 2. 2| fixed pressure 
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A ( + A r = 


§p.2.2 fixed flow rate 
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= 0, 
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§p.2.2 fixed temperature 


6 


TO 


= 0.01 


7T = 


0.99 


A r = 


§p.2.3| nonzero terminated velocity; 














"leaky flow" 



The condition A^ = means that the density does not oscillate at the lower boundary. The 
condition A^ + A T = stands for constancy of the mass flow rate pv at the lower boundary. The 
physical interpretation of these two choices has a subtle complication: in the present treatment 
the density grows indefinitely as £ —>■ 0. However this problem would be alleviated in a modified 
geometry (e.g. taking account of lateral flow from the sides of the accretion column) which will 
be investigated in future studies. (The present model applies to white dwarf accretion shocks so 
long as lateral spreading of accreted gas is much slower than infall at the shock, and buildup of 
material at the base is slow on the cooling and oscillatory timescales.) 

Other physical lower boundary conditions are related to the thermal structure of the post- 
shock flow. The electron and ion temperatures must be equal at the white dwarf surface, which 
implies that A c — A„- = is a necessary boundary condition. However this condition is satisfied 
for every choice of 5, since the functional form of the electron-ion energy exchange F ultimately 
exceeds the bremsstrahlung cooling Ab r in the high-density, low-temperature region near the lower 
boundary (see e.g. Saxton & Wu 2001). Another choice would be to fix the flow temperature at 
the lower boundary, assuming that the gas at the base of the flow is in thermal equilibrium with 
the white dwarf, and that the white dwarf surface is a perfect heat-bath/heat-sink. Constancy 
of T oc P/p implies a condition of A^ — A^ = 0. The lower-boundary temperature is zero in the 
stationary solution for the present example, so this choice of condition on the perturbed variables 
describes the lower boundary as being a perfect node in the oscillation of the temperature profile. 

3.2.3 nonzero terminated velocity 

This paper also considers an example in which the lower boundary has a non-zero flow velocity 
relative to the shock. This condition describes a system where cooling at the base does not 
compress the gas indefinitely, and the the accretion shock moves outwards from the white dwarf 
surface (until non- linear effects become important or the shock becomes indistinct). Alternatively 
the boundary condition could describe a surface at which post-shock material is removed from 
the system at a nonzero rate, like flow into a "leaky" bucket. The illustrative choice considered 
here is To — > 0.01 rather than zero, and A r = at the lower boundary. In other words, the shock 
propagates away from the white dwarf surface at a velocity that is 0.01 times the free-fall velocity 
of pre-shock matter, or else gas is removed from the accretion stream where it contacts the stellar 
surface at incident velocity O.OIvq. 
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4 Stability properties 



4.1 eigenf unctions in the standard case 

In the "perfect stationary wall" case, there exist amplitude nodes and antinodes in the eigen- 
functions for perturbed density (|A^|), and also for electron pressure (|A e |) in cases when the 
two-temperature effects are important. The number of nodes depends on the harmonic number of 
the mode. Phases of the eigenfunctions change abruptly at the node locations, and wind gradually 
in intervening regions. 

Apart from the perturbed velocity (A T = by the boundary condition) the perturbed variables 
are not explicitly determined at the lower boundary, but have finite amplitudes within an order 
of magnitude of their typical amplitudes throughout the whole post-shock region. The perturbed 
electron pressure and total pressure variables always meet at the same value at the lower boundary, 
regardless of the boundary condition, because the electron and ion temperatures equilibrate at 
the lower boundary. For comprehensive discussion of the eigenfunctions, see Saxton (1999) and 
Saxton & Wu (2000). 

4.2 sensitivity of eigenfunctions 

Saxton (1999) §5.4.1-5.4.2 considered the relationships between eigenvalues and local features of 
the corresponding eigenfunctions. The relevant findings were as follows. 

For arbitrary alterations of the lower boundary condition resulting in a small offset of the 
eigenvalue, \A8\ = 0.01, the eigenfunctions alter significantly in the region ^0.05 near the lower 
boundary, but remain essentially unchanged throughout the majority of the post-shock region 
(0.05^,^ < 1). In the sensitive base region, offsets from the "perfect stationary wall" condition 
yield steeper gradients in the amplitudes and phases of all the A- variables. The amplitudes |A^| 
and |A T | may differ their ordinary values by over an order of magnitude; the pressure amplitudes 
| j and | A c differ by factors of a few. 

At more extremely different boundary conditions, resulting in eigenvalues that are intermediate 
between modes of the "perfect stationary-wall" case, the sensitive, base region (£,$,0.03) with 
steep gradients in the A-functions persists. Amplitudes are enhanced in regions up to £^0.85, but 
this is not as significant as the efffects in the base region. The number of antinodes and nodes 
resmbleare comparable to those of the next higher modes of the "perfect stationary-wall" case, 
but their positions are completely different. 

Thus small offsets of the global frequency and stability properties alter the local behaviour of 
the lowest few percent of the post-shock flow, but have no significant effects in regions of the flow 
further from the lower boundary. 

4.3 eigenvalues: stability and frequencies 

The previous subsection reviewed how eigenfunctions of the hydrodynamic variables generally 
respond to alterations of lower boundary conditions resulting in small and large alterations of the 
frequencies and stability properties. This section details the effects on the eigenvalues, the mode 
frequencies and stabilitiy properties, of specific alternatives for the lower boundary condition. 

4.3.1 "perfect stationary-wall" 

The general stability and frequency results for the standard "perfect stationary- wall" lower bound- 
ary condition are as follows (see e.g. Chevalier & Imamura 1982, Imamura et al. 1996, Saxton et 
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al. 1997, 1998 and Saxton & Wu 1999, 2000). When bremsstrahlung emission is the only cooling 
process the fundamental mode is stable, and the overtones are unstable, with instability 5 R tend- 
ing to increase with harmonic number. When two-temperature effects are unimportant (ip e i > 1), 
increasing the cyclotron cooling efficiency (e s ) stabilises each mode. Some modes stabilise more 
readily than others and the trend of 6 R increasing with harmonic number breaks down more ex- 
tensively as cyclotron cooling becomes more efficient. In the absence of two-temperature effects, 
the oscillatory frequencies are quantised like those of a pipe open at one end (Toth & Draine 1993; 
Saxton et al. 1997, 1999). Increasing e s reduces the frequency intervals of between consecutive 
modes (Saxton et al. 1997; Saxton & Wu 1999). When two-temperature effects are important 
(ipd^l), greater e s may actually destabilise modes (Imamura et al. 1996) and the frequency se- 
quence becomes more like that of a doubly-open or doubly-closed pipe (Saxton & Wu 1999). 

Small offsets of the 5 values may produce large changes in the A-profiles in the region near 
the lower boundary. The affected region is only a few percent of the entire post-shock flow. For 
larger 5 offsets the A-profiles are qualitatively similar to the eigenfunctions of the nearest higher 
mode, except that: (1) the perturbed variables change abruptly near the lower boundary, as for 
small offsets; (2) throughout most of the post-shock region there is a general increase of amplitude 
with decreasing £; and (3) the nodes and antinodes are located at different positions. Thus 
systems with eigenfrequencies and stability properties very different from the stationary-wall case 
can nonetheless have similar local oscillatory responses throughout most of the post-shock region 
except for a small base region. 



4.3.2 modified stationary wall 



Table ^ presents eigenvalues for the alternative boundary conditions described in § |3.2.2| : constant 
pressure, flow rate, density or temperature. Compared to the "perfect stationary-wall" case, all of 
the modified stationary- wall conditions yield complex 5 values that differ by only a few percent. 
The changes of 5 are smaller than or comparable to the idealised "small offset" cases (§ [4.2| above; 
also Saxton 1999, §5.4.1)). 

Thus the alteration of the lower boundary condition on perturbed variables has no significant 
effect on the global stability properties and oscillatory frequencies of the modes. The nature 
and efficiencies of the radiative cooling processes and the electron-ion energy exchange have a 
greater influence on the frequencies and stability properties of the oscillating shocks than the lower 
boundary conditions. Observed frequency ratios and damping characteristics of shock oscillations 
(e.g. in accretion onto magnetic white dwarfs) is expected to be practically independent of the 
lower boundary condition. (However the lower boundary may still have significant effects upon 
the modulation of bremsstrahlung and cyclotron luminosities, as discussed in §5.) 



4.3.3 nonzero terminated velocity 

If the flow velocity terminates with a nonzero value, tq = 0.01, at the lower boundary (§ [3.2.3|) , 
the oscillatory frequencies are essentially unchanged, and the stability properties change slightly. 
For the cases studied, when bremsstrahlung emission is the only cooling process the difference 
between zero-velocity termination and nonzero-velocity termination is a few times 0.001 in both 
5 R and <5 t ; for e s = 100 the chang few times 0.01. Changing from zero-velocity termination 

to nonzero-velocity termination at the lower boundary made every mode more stable. Higher- 
order modes are more significantly affected, and the effects are greatest when e s is higher. As 
seen in Figure |I|, the sequence of the modes' stability properties is qualitatively different for the 
zero- and nonzero-terminated velocity boundary conditions in the cyclotron-dominated extreme 
K Vei, e s ) = (0.5,0.5,100). 
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Table 2: Eigenvalues for four alternative choices of the lower boundary condition in the perturbed 
variables, but with the standard conditions on the stationary solution (r = 0, 7r = 1, n c = \ at 
£ = 0). Values shown are for the first six harmonics, n — 1 . . . 6. 
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The increased stability of modes in the cases of cyclotron-dominated cooling may occur because 
the conditions at the lower boundary are not as dense as in the "perfect stationary-wall" case. 
Since the highest densities near the lower boundary are less extreme, the bremsstrahlung cooling 
dominates over a smaller region at the base of the post-shock flow. As bremsstrahloung radiative 
cooling is responsible for the thermal instability that drives the shock oscillations, it is qualitatively 
sensible that the "nonzero terminated velocity" case has more stable modes than the "perfect 
stationary-wall" case. 

5 Luminosity response 

Oscillations of the post-shock flow structure modulate the emitted radiation. Rapid ~ 1Hz shock 
oscillations were observataionally inferred in accretion shocks of white dwarfs in some AM Her- 
culis (AM Her) systems by optical-infrared observations of the oscillating cyclotron luminosities 
(Middleditch 1982; Imamura & Steiman-Cameron 1986; Larsson 1987, 1989; Ramseyer et al. 1993; 
Middleditch, Imamura & Steiman-Cameron 1997). However the shock oscillations of these sys- 
tems may not produce comparable oscillations of bremsstrahlung luminosity. For some systems 
the X-ray luminosity responses are subject to upper limits of a few percent (Wolff et al. 1999; 
Beardmore & Osborne 1997). It is therefore interesting to predict the oscillatory responses of the 
bremsstrahlung and cyclotron total luminosities implied by the analytically obtained modes, and 
their phasing properties (Saxton & Wu 2001). 

Wolff et al. (1991) and Wood et al. (1992) carried out numerical investigations of noise-drive 
QPOs in cases where the magnetic field strength is low enough that the cyclotron radiative cool- 
ing does not dominate. They found that modes that are suppressed in linear analyses may be 
sustained by accretion noise and show significant luminosity responses. The luminosity responses 
were found to be linear, and bremsstrahlung and cyclotron amplitudes were comparable although 
bremsstrahlung signatures of the modes are broader and therefore less distinct in the frequency 
power spectrum. 

The present work extends the previous analysis to consider how the luminosity oscillations 
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may depend on the character of the lower boundary. The relative suppression ofrapid oscillations 
in the X-rays may probe the nature of the boundary at the bottom of the post-shock accretion 
flow onto the white dwarf in AM Her systems. 

The different radiative cooling processes and the electron-ion energy exchange are strong in dif- 
ferent regions of the post-shock flow. Bremsstrahlung cooling depends on a high power of density, 
and therefore it peaks in the cold, high-density region near the lower boundary. Cyclotron cooling 
depends on a high power of temperature, and therefore most of the cyclotron emission comes 
from the hot region near the shock. Therefore the local oscillations of hydrodynamic variables 
in the region near the lower boundary are effective at modulating the bremsstrahlung luminosity, 
whereas the oscillatory behaviour of the flow near the shock is most effective at modulating the 
cyclotron luminosity. Thus the bremsstrahlung luminosity and the cyclotron luminosity have dif- 
ferent responses to the shock oscillation of any given mode. Furthermore, because different modes 
have different profiles of nodes, antinodes and phases between the upper and lower boundaries, 
the oscillatory responses of the bremsstrahlung and cyclotron luminosities must differ between 
modes. 

As shown in Saxton (1999), the oscillatory eigenfunctions are sensitive to changes of the lower 
boundary condition in the region close to the boundary (£,$,0.03) but insensitive in the rest of the 
post-shock flow (0.03^,^ < 1). Since bremsstrahlung emission peaks near the lower boundary and 
cyclotron emission peaks near the shock, it is to be expected that the integrated bremsstrahlung 
luminosity is sensitive to alterations of the lower boundary condition, but the cyclotron luminosity 
is not. The bremsstrahlung luminosity response is additionally sensitive to the lower boundary 
condition by way of the large amplitudes in the density eigenfunction A^ that can occur at low 
£ for some S frequency values. The calculated luminosity responses shown in Table |3| confirm 
this prediction when the "perfect stationary-wall" condition (A r = 0) at the lower boundary 
is replaced by the alternative conditions described in § ^.2.2| and § |3.2.3j . The bremsstrahlung 
luminosity amplitudes and phasings change significantly under alterations of the lower boundary 
condition, but the cyclotron luminosity response is essentially unaffected. 

Because of the insensitivity of cyclotron luminosity responses to the lower boundary condition, 
the following discussion concentrates on the bremsstrahlung response alone. 

5.1 "perfect stationary- wall" 

As discussed for the choice of a "perfect stationary- wall" lower boundary in Saxton & Wu (1999), 
the modes' luminosity responses were found to be insensitive to (a s ,ip e i) when bremsstrahlung 
cooling dominates (e s = 0). However there are no obvious relationships between the luminosity 
responses and stability properties of the modes for diverse cases of the system parameters. It 
was found that the cyclotron luminosity oscillation tends to lag the bremsstrahlung luminosity 
oscillation by approximately 0.67T radians for the fundamental mode, but the overtones lack obvious 
relationships between their bremsstrahlung and cyclotron responses. 

5.2 fixed-density or fixed-temperature 

The bremssttrahlung luminosity responses are identical to within a few percent for the cases 
of fixed-density (Af = 0) and fixed-temperature A^ — A^ = 0) boundary conditions. Thus in 
this sense an absence of oscillation in density at the lower boundary is almost equivalent to an 
absence of oscillation in temperature. For each choice of (<r s , ip e i, e s ) there is an approximately 
constant ratio between the |L bril |/eL br o values in the fixed-density and fixed-temperature cases 
as compared to corresponding modes in the "perfect stationary wall" case (see Figure 0). For 
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Figure 2: Comparison of bremstrahlung luminosity oscillatory amplitude for cases with "perfect 
stationary- wall" (A T = 0, solid lines) and fixed-density (A^ = 0, dotted lines) lower boundary 
conditions. Values for the first six harmonics are marked with stars, squares and diamonds for 
e s = 0, 1, 100 respectively. The other system parameters are (a s , ip c \) = (0.5, 0.5). For given system 
parameters, the |Lb ri i|/£Lbr,o values in the fixed-density case are approximately proportional to 
corresponding values in the "perfect stationary-wall" case. 
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(<t s , ipei, e s ) = (0.5,0.5) the "perfect stationary-wall" and fixed-density/temperature |Lbr,i| values 
are in the approximate ratio of 1.0 for e s = 0, 1.0 — 1.2 for e s = 1, and 1.3 — 1.7 for e s = 100. 
The phases of the bremsstrahlung luminosity responses are approximately equal for each mode in 
the fixed-density/fixed-temperature cases and the "perfect stationary- wall" case. For the studied 
modes where (cr s ,-?/>ei) = (0.5,0.5), the phase differences are ^0.027r radians for e s = 0, ;$0.037r 
radians for e s = 1, and ^,0.047r radians for e s = 100, 

5.3 fixed-pressure 

With the present choices of system parameters, the bremsstrahlung luminosity of each mode is 
greater under the fixed-pressure condition (A^ = 0) than under the "perfect stationary- wall" condi- 
tion (A T = 0). Modes of higher harmonic number tend to have higher amplitudes of bremsstrahlung 
luminosity response. (One exception is the n = 3 mode when (<r s ,^ e i,e s ) = (0.5,0.5, 100), which 
is also unusual because the instability eigenvalue 5 K departs from the usual trend of 5 R increas- 
ing with harmonic number n.) The fixed-pressure and "perfect stationary- wall" lower boundary 
conditions also produce different bremsstrahlung luminosity phases. As shown in the top-left and 
bottom-left blocks of Table [|, the bremsstrahlung luminosity phases for the first six harmonics 
differ by ~ 0.37T — 0.87T for (o" s ,^ ei ) = (0.5,0.5) and e s = 1, 100. The change of boundary condi- 
tion makes less difference (^,0.37r in phase) for the same modes when cyclotron cooling is absent 
(C = 0). 

5.4 fixed flow-rate 

In the cases with fixed flow-rate boundary condition (A^ + A T = 0) the bremsstrahlung luminosity 
response amplitudes generally decrease with harmonic number, n. The exceptions, for the choices 
of system parameters studied here, are the modes which depart from the trend of <5 R increasing 
with harmonic number, n. There is no clear relationship between the phases of bremsstrahlung 
luminosity oscillation in the "perfect stationary-wall" and fixed-flow cases, however the phase 
of each mode in the fixed-flow case is approximately in antiphase to that in the corresponding 
fixed-pressure case. 

5.5 nonzero terminated velocity 

When cyclotron cooling is absent (<r s , ip ei , e s ) = (0.5,0.5,0), the bremsstrahlung luminosity re- 
sponses for the nonzero terminated velocity condition (r = 0.01, A r = 0) are almost the same 
as for the "perfect stationary-wall" condition (r = 0, A T = 0). The corresponding phases are 
approximately equal (to within ^0.0057r radians, see Table |^, bottom row). The amplitudes agree 
to within a few percent. 

When bremsstrahlung and cyclotron cooling are equally efficient at the shock (cr s , e s ) = 
(0.5,0.5,1), zero- and nonzero- velocity termination have different luminosity responses for the 
lowest three modes, but approximately equal luminosity responses for the higher-order modes. 
When cyclotron cooling dominates, (cr s , ipei, e s ) = (0.5,0.5,100), the velocity termination affects 
the luminosity responses for all the modes studied. 

It is not clear why nonzero velocity termination at the lower boundary affects luminosity re- 
sponses of low-harmonic modes more than higher modes. Higher modes have Ab r profiles with 
more "nodes" and "antinodes" , which should be expected to make them more sensitive to altered 
conditions at many locations in the post-shock flow; lower modes have fewer amplitude extrema 
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and lesser spatial gradients of phase. It can be speculated that for modes of higher harmonic num- 
ber, the integration in ([14]) involves a proportionately smaller contribution in the low-£, compared 
to the case of lower modes. 

5.6 Comparison with numerical results 

In numerical studies of noise-driven shock oscillations, in low-field (low e s ) cases, Wood et al. (1992) 
determined the lag of cyclotron luminosity response relative to the bremsstrahlung luminosity 
response, which is comparable to the difference $ cy — $br in the notation of the present paper. For 
cases with comparable equilibrium luminosities of cyclotron and bremsstrahlung (L cyj0 /£br,o = 0.5) 
the cyclotron phase lag was — 0.87T radians for the fundamental (n = 1) mode, and — 0.2it for the 
first overtone (n = 2). This result is approximately consistent with the typical lag ~ — 0.67T for 
the fundamental mode in the fixed-density, fixed-temperature and fixed- velocity lower boudnary 
conditions for cases with e s = 1 in Table ^ but inconsistent with fixed-flow-rate and fixed-pressure 
conditions. Comparison of the cyclotron luminosity lags for the n = 2 mode may not be meaningful 
because in the analytic results, the n = 2 mode's phasing is sensitive to the system parameters 
as well as the lower boundary condition. Wood et al. (1992) also caution that phasing may differ 
between free and driven oscillations; by assumption, our analytically represented oscillations are 
physically free. 

As Imamura (1985) recognised, the discrepent results from numerical investigations of radiative 
shock instabilities in accreting white dwarfs may be largely affected by the numerical treatment of 
the cold, high-density, gas near the lower boundary, where rapid cooling causes several of the fluid 
variables to experience sharp gradients. In the Lagrangian numerical hydrodynamic calculations 
of Imamura, Wolff & Durisen (1984) the grid zones accumulated at the lower boundary were 
eliminated by coalescing the cells. In Imamura (1985) the cluttered grid zones were instead 
eliminated after cooling below a threshold temperature. For the latter method the shock stabilities 
were less sensitive to the procedural details than in the previous method. It is possible that 
methods involving mergers or eliminations of material accreted at the base of the flow have an 
effectively open lower boundary condition. If the merger /deletion process still conserves mass, 
energy and momentum then it still may effectively delete structural information such as the 
ideally steep gradients of velocity and density near the lower boundary. The importance of those 
gradients is shown by the different oscillatory responses of the luminosities, between the "perfect 
stationary- wall" and "nonzero terminated velocity" boundary conditions. 

The Eulerian calculations (e.g. Langer, Chanmugam & Shaviv 1981, 1982; Chanmugam, Langer 
& Shaviv 1985) necessarily introduced a spatial filter to remove noise and numerical instabilities. 
This process may also affect the thermal instability by smoothing away oscillatory structures, such 
as those in the sensitive region immediately above the lower boundary, that vary on finer spatial 
scales than the filtering resolution (e.g. nodes and phase jumps like those in the eigenfunctions of 
higher modes in Saxton & Wu 2001). 

6 Conclusions 

The oscillatory and stability properties are not significantly affected for modifications of the 
stationary-wall boundary condition involving the perturbed hydrodynamic variables: for con- 
stant pressure, flow rate, density or temperature at the lower boundary. The cases considered 
here have eigenvalues which are effectively identical to the "perfect stationary- wall" results, or at 
most the "small offset" profiles which only affect the A-functions is regions very close to the lower 
boundary. 
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Table 3: Bremsstrahlung and cyclotron luminosity responses, for the first six modes, for various 
choices of the system parameters and lower boundary condition. For the blocks of values from left 
to right, top to bottom, the respective lower boundary conditions is A,r = 0, A^ + A T = 0, A^ = 0, 
A^ — A^ = 0, (A T = 0, r = 0) and (A T = 0, r = 0.01). The amplitudes are normalised by e and 
the phases are multiples of n radians. 
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If the stationary solution terminates at a lower boundary with nonzero velocity, the oscilla- 
tory frequencies are essentially unchanged from the "perfect stationary-wall" case. However the 
nonzero-terminated velocity modes are slightly more stable than the corresponding modes of the 
"perfect stationary-wall," and the higher modes are more greatly affected. 

The cyclotron luminosity response is independent of the lower boundary condition as this 
emission is peaked near the shock, where jump conditions specify the hydrodynamic variables and 
the lower boundary conditions have little influence. The bremsstrahlung luminosity response is 
a sensitive probe of the conditions at the lower boundary, because bremsstrahlung emission is 
peaked near that boundary and the perturbed hydrodynamic variables in this region are locally 
very sensitive to changes in the lower boundary conditions. However the fixed-density and fixed- 
temperature lower boundary conditions have modes with identical luminosity responses; the reason 
for this equivalence is unclear. 

Termination of the stationary solution at nonzero velocity produces effectively identical lumi- 
nosity responses as the same modes with zero velocity termination, for all modes when cyclotron 
cooling is unimportant (i.e. when e s is small). As e s increases, the velocity termination makes 
a difference to the luminosity responses of the lowest-order modes but not higher-order modes. 
When e s is sufficiently large, nonzero velocity termination affects the luminosity responses of all 
modes studied. For modest two-temperature effects, the bremsstrahlung luminosity responses of 
higher-order modes tend (with few exceptions) to be greater under the fixed-pressure lower bound- 
ary condition, and lower under the fixed flow-rate lower boundary condition. However there is no 
universal pattern between 5 and the luminosity responses for general (a s ,ip e i). For a given mode 
and given system parameters, the bremsstrahlung responses in the fixed-pressure and fixed-flow 
cases are approximately in antiphase. 

Past numerical simulations reproduce some of what may be the important qualitative features 
of these thermally unstable radiative shocks, though it is not clear which, if any, of the numerical 
treatments provides superior predictions about the stability properties and frequencies of the 
shock oscillations. All the treatments to date may lack sufficient spatial and temporal resolution 
in the lower boundary region, where the oscillatory instability may originate. Indeed treatments 
may alter the fundamental character of the boundary, (e.g. the removal of grid zones or accreted 
material may provide a more open boundary than would be realistic for the fixed wall of the white 
dwarf surface). 
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